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Abstract 

We use accurate quantum mechanical calculations to analyze the effects of parallel electric and 
magnetic fields on collision dynamics of OH(^n) molecules. It is demonstrated that spin relaxation 
in ^He-OH collisions at temperatures below 0.01 K can be effectively suppressed by moderate 
electric fields of order 10 kV/cm. We show that electric fields can be used to manipulate Feshbach 
resonances in collisions of cold molecules. Our results can be verified in experiments with OH 
molecules in Stark decelerated molecular beams and electromagnetic traps. 
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I. INTRODUCTION 



The rapid progress in the research field of cold molecules holds greatpromise for new 

and important discoveries at the interface of physics and chemistry y,iiQ,y,y,3. The 

remarkable properties of cold molecules allow for their application in quantum information 
processing, condensed-matter physics, physical chemistry, and precision spectroscopy. While 
some of these applications rely on specific properties of molecular ensembles, others exploit 
the diversity of molecular electronic states and energy level structures. For example, re- 
cently proposed schemes for quantum computing with polar molecules in optical lattices 
make use of intermolecular dipole-dipole interactions and molecule-field couplings to enable 

n 

communication between quantum bits pj. The experiments aiming to measure the electric 
dipole moment of the electron based on the spectroscopic study of cold YbF radicals exploit 
the relativistic distortion of molecular orbitals [2i]. The fine and hyperfine perturbations in 
the spectra of ^11 molecules allow for the development of new frequency standards and tests 



of physics beyond the Standard Model [3,|4] 

Particularly noteworthy are the experiments probing inelastic energy transfer and chem- 
ical reactions of cold molecules 5|, |6|, |7|. The measurements of molecular collision properties 
at low temperatures are important for the development of buffer gas cooling techniques, 
which rely on coUisional thermalization of molecules in a cell filled with cryogenic ^He gas. 
This technique produces molecules at temperatures between 0.1 K and 0.5 K and allows 
for further evaporative cooling after removal of the buffer gas. In order to permit obser- 
vations and evaporative cooling, molecules must remain trapped for a long time, which is 



possible if the ratio of elastic to inelastic spin-changing collisions exceeds 10^ 
experimental [ol, |lo| and theoretical 



Previous 



ll| work has shown that this requirement is generally 



satisfied for S-state diatomic molecules with large rotational constants. However, it is un- 
clear whether molecules in electronic states other than S may be sympathetically cooled in 
a buffer gas cell. The mechanism of spin relaxation in collisions of molecules in ^11 electronic 
state remains unknown. 

The OH radical is one of the first molecular species that were cooled and trapped at 
milliKelvin temperatures 



12 



The ground electronic state of OH is of ^H 
symmetry, and the energy levels of the molecule depend linearly on the strength of an applied 
electric field above ~5 kV/cm. As a result, the OH radicals can be efficiently decelerated. 
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trapped, and manipulatec 



the laboratory [5|, Il2| 
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using moderate static or time-varying electric fields available in 



14| . A variety of novel experiments with trapped OH molecules 



has been reported. The ability to fine tune the collision energy of a Stark-decelerated beam 
allowed Meijer and co-workers to study threshold behavior of Xe-OH rotationally inelastic 
scattering at collision energies as low as 50 cm~^ [5]. High-precision spectroscopy of cold 
trapped OH was used to study the time evolution of the fine structure constant [4i], measure 
the lifetimes of vibrationally excited states 1^ , and determine the rates of optical pumping 
due to blackbody radiation . Sawyer et al. have recently reported measurements of cross 
sections for elastic and inelastic collisions of magnetically trapped OH molecules with He 
atoms and D2 molecules at kinetic energies of 60 cm^^ and above [6|. 

Theoretical studies of low-energy collisions of OH molecules have been reported by several 



groups 
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171 . Il8l . Il9|. Avdeenkov and Bohn studied ultracold collisions of OH 



molecules [15|, ll6|] and discovered weakly-bound dimer states supported by the dipole-dipole 



interaction forces in the presence of an external electric field [ifii]. Ticknor and Bohn found 
a significant suppression of inelastic relaxation rates for the low-field-seeking states of OH in 
a magnetic field {itI . Lara et al. analyzed the effects of non-adiabatic and hyperfine effects 
on field-free Rb-OH collisions [18]. Their results indicated that sympathetic cooling of OH 
molecules by collisions with Rb atoms might be challenging due to large inelastic loss rates. 



Gonzalez-Sanchez, Bodo, and Gianturco 



19j considered field-free collisions of rotationally 



excited OH molecules with He atoms and found sharp propensity rules for rotational and 
A-doublet changing transitions at ultracold temperatures. 

Here, we present a theoretical analysis of OH(^H) collision dynamics in combined electric 
and magnetic fields. We have previously demonstrated that spin-changing collisions of S- 
state molecules can be efficiently manipulated by superimposed electric and magnetic fields 
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15 
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. Building on our previous work [20|, l21j and the results of Bohn and co-workers 
171], we develop a rigorous quantum theory of collisions between ^H molecules and 
structureless atoms in external fields and calculate the dependence of the cross sections 
for He-OH collisions on electric and magnetic fields. Our results suggest that collisions of 
OH molecules with He atoms can be efficiently manipulated with the external fields. In 
particular, we demonstrate an efficient mechanism for suppression of spin relaxation in ^H 
molecules with electric fields. 



3 



II. THEORY 



The quantum mechanical formahsm for colhsions of diatomic molecules in ^11 electronic 
states in the absence of external fields has been presented by several authors (see, e.g., Refs. 



19|, |23|, |2J] ) . Here, we focus on the theoretical aspects relevant for incorporating the effects 



of electromagnetic fields in scattering calculations. Section IIA presents the discussion of the 
influence of the electric and magnetic fields on the energy level structure of ^11 molecules. 
In Sec. IIB, we discuss the Hamiltonian of the collision complex and the coupled- channel 
representation of the scattering wave function. A derivation of the matrix elements for the 
interaction potential operator in the basis of scattering states is presented in Sec. IIC. 



A. The OH molecule in superimposed electric and magnetic fields 



The Hamiltonian for a H molecule such as OH can be written as [25|, l26|, |2j 

-f^mol = -f^rot + -f^so + Ha + He + Hb, (1) 
where i^rot is the angular part of the rotational kinetic energy 28 1 



Hrot = i?e(i'" - - S'"")', (2) 

Be is the rotational constant, J^^ = N^'^ + L^^ + S^^ is the total angular momentum, 
iV^^ is the rotational angular momentum of the nuclei, L^^ is the electronic orbital angular 
momentum, and S^^ is the electron spin. In Eq. ([1]), we have neglected the hyperfine 
interaction due to the nuclear spin of H. The hyperfine interaction constant of OH is an 
order of magnitude smaller than the A-doublet splitting, and the hyperfine effects may alter 
collision dynamics at temperatures below 4 mK [17]. The angular momentum operators 
in Eq. ([2]) are defined in the space-fixed frame. However, the symmetry properties of the 
electronic wave functions are most conveniently exploited in the molecule-fixed frame, with 
the z-axis oriented along the OH bond. The row vector of molecule-fixed angular momentum 
operators can be defined as J^^ = J^^R{a, (3,0), where R is the matrix of direction cosines 
and (a,/?) are the Euler angles which specify the orientation of the diatomic molecule in 
the space-fixed coordinate system. Alternatively, one can define the column vector J^^ = 



R{a,p,0)J 



MF 



28j. The molecule-fixed angular momentum operators do not commute, and 



the choice of convention affects the products of operators that arise upon transforming 



Eq. ([2]) to the molecule-fixed frame. It is easy to show that the frame transformation 
adopted in this work leads to J^^ ■ S^^ = ■ S^^, whereas that of Ref. 3] leads to 
jSF . gSF _ gUF . jMF r-|nj_^g matrix elements of H,.ot are independent of the convention. In 
the following, we will omit the superscript "MF" . 

The second term in Eq. ([T]) is the spin-orbit (SO) interaction 



^so = AL-S, (3) 

where A is the SO interaction constant. The remaining terms in Eq. ([T]) account for the 
effects of A-doubling and the interactions with static electric and magnetic fields (explicit 
expressions for -^e, and Hb are given below). The energy levels of a ^11 molecule can 
be evaluated by diagonalizing the Hamiltonian ([1]) in Hund's ) basis 

/2/+i\V2 

\jMn)m = i^-^j Diua,mm, (4) 

where M and Q are the projections of J onto the space-fixed and molecule-fixed quantization 
axes, -D^/j^ is the Wigner D-function, and |AS) is the electronic wave function. The molecule- 
fixed projections of L and S are denoted as A and S. For a ^11 electronic state, they take 
the values A = ±1 and S = ±i. The off-diagonal matrix elements of the SO interaction 
between the ^11 state and the excited electronic states give rise to the A-doubling effect 



261 ] described below. After neglecting the cross terms, Eq. ([3]) may be written as 

^so = AL^S,, (5) 
where the subscript z refers to the molecule-fixed projections of the angular momentum 



operators. The A-doubling is described by the effective Hamiltonian [25|, [26 1 



= le-'^^'l-qJl + {P + 2q)jJ+] + \e^'^[-qJl + (p + 2g) J_^_], (6) 

where J± = Jx T iJy and S± = Sx ± iSy are the ladder operators, is the azimuthal 
angle of the electron in the molecule-fixed frame, and p and q are the phenomenological 
A-doubling parameters. Using the phase convention for the electronic wave functions, (A = 
±l|e^^*'^|A' = =1=1) = —1, the matrix elements of the Hamiltonian ([6]) can be written as 

(A|^a|A') = i {SA',-i5A,i[qJl -ip + 2q)J+S+] + (5AM(5A,-i[gJ! - (p + 2g) J_^_]) (7) 
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The interaction with the magnetic field of strength B has the form 

Hb = fioBii + 25) ■ 5, (8) 

where B is the unit vector in the direction of the external magnetic field. To first order, the 
interaction of the molecule with the dc electric field can be written as 

He = —Ed cos X, (9) 

where x is the polar angle of the molecule in the space-fixed frame, E is the electric field 
strength, and d is the permanent electric dipole moment of the molecule. Here, we assume 
that both the electric and magnetic fields are oriented along the space-fixed 2;-axis. The 



29|. 



more general case of crossed electric and magnetic fields is considered elsewhere 
It is convenient to use parity-adapted Hund's ) basis functions 

IJMQe) = ^{\jMn)\A = = n - 1) + e{-y-^/^\jM-n)\A = = -0 + 1)} , 

(10) 

where Q = and the parity index e = ±1 characterizes the inversion symmetry of the 
basis functions [e{—)'^~^^'^ = 1 for e-parity states and —1 for /-parity states]. Note that in 
the parity-adapted basis, > and the quantum number A does not have a definite value. 
Expanding the rotational kinetic energy in terms of the ladder operators and using Eq. (I5l), 
we obtain 

Hrot + Hso = B,[J^ - 2Jl - J+5_ - J_5+ + {A + 2B,)lJ,], (11) 

where we have omitted the terms and 5^ which would only result in an overall energy shift. 
The matrix elements of the rotational and spin-orbit Hamiltonians can now be evaluated in 
the parity-adapted basis ( fTOl) . They have the form 

{JMne\Hrot + Hso\J'M'n'e') = 6,,>6jj,6MM'{Be[J{J+l)-2n^]6m' + {A + 2B,)^^ 

- B,[6^^n,^,a^{J',Q')a4S,n' - 1) - 5n,n'+i«+(^', ^O^+l-^, 0' - 1)]}, (12) 

where a±(J, 0) = a/J(J + 1) - 0(0 ± 1). Combining Eqs. ([7]) and (flO]), we obtain the 
following compact expression for the A-doubling matrix elements 

{JMQe\Hj,\J'M'Q'e') = ^<5,,.5jj,We(-)''"'/'[ga-( J', n')a-( J', H' - l)%,2-f^, 

- (p + 2g)a_( J', n')a+{S, 0' - l)5n,i-n']- 

(13) 
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In order to evaluate the matrix elements of the Zeeman Hamiltonian in the basis (fTOl) . it is 
necessary to transform the operator ([8]) to the molecule-fixed frame 

1 

Hb = fioBY,iLg + 2SM*{a,P,0), (14) 

g=-l 

where only the q = molecule-fixed component of L survives on the right-hand side. Eval- 
uating the integrals over the product of three Wigner D-functions, we find 



{JMQe\HB\J'M'Q'e') = /xq^^^'W' (-)''''"''' [(2 J + 1)(2J' + 1)]^/^ 



J 1 J' 
M -M' 



X 



r- -, I J ^ J' \ r- -, J 1 J' 

V2a+iS,n' -1) \ _ _ - V2«_(S,f]'- 1) 

n -1 -n' J \ n 1 -n' 

+(2n-i; 



J 1 J' 
n -n' 



(15) 



The matrix elements of the interaction with electric fields have a similar form 
{JMne\HE\J'M'n'e') = -Ed5,,.,,5MM'i-)^'~^'[{2J + 1)(2J' + 1)]^/^ 

x| ' \ (16) 

M -M' j \ -fi' 

This expression shows that electric fields couple the states of the opposite inversion parity. 



B. Collision dynamics 

The He-OH(^n) collision complex can be described by the Jacobi vectors R - the separa- 
tion of He from the center of mass of OH and r - the internuclear distance in OH. The angle 
between the vectors is denoted by 9. In the following, it will be convenient to use the unit 
vectors R = R/R and f = r/r, where R = \R\, f = \r\. The Hamiltonian of the collision 
complex can be written in atomic units as [ij 
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23] 



where fi is the reduced mass of the ^He-OH system, i is the orbital angular momentum for 
the collision, V{R, r, 6) is the electrostatic interaction potential, and Hmoi is the Hamiltonian 



of the OH molecule in the presence of external electric and magnetic fields (see Sec. IIA). 
We assume that the internuclear distance of OH is fixed at the equilibrium value of 1.226 A. 
The wave function of the He-OH collision complex \1/ satisfies the Schrodinger equation at 
a total energy E, and can be expanded over the complete coupled-channel basis 

'^ = ^Y.^P(^^MR^r), (18) 

where ippiR, r) are the angular basis functions. The fully uncoupled angular basis set can be 
defined as a direct product of the parity-unadapted Hund's ) functions and spherical 

harmonics 

\JMn)\^.T)\^mi), (19) 

where the spherical harmonics \l'm^) = Yf„n{r) describe the orbital motion of the He atom 
around the OH fragment. To be consistent with spectroscopic nomenclature, it is convenient 
to use a slightly modified basis given by 

\JMCLe)\lmt), (20) 

where |JiVff2e) are Hund's ) basis functions of definite parity (fTOj) . The basis sets 

([T9|) and (|20|) are related by a unitary transformation, and are equivalent. 

Substituting the coupled-channel expansion (fT8|) into the Schrodinger equation, we obtain 
a system of coupled second-order differential equations 

Fp{R) = 2^iY,{MR,f)\V{R,e) + ^ + H^^,\ijp,{R,f))Fp,{R). (21) 

In order to solve these equations, it is necessary to evaluate the matrix elements on the right- 
hand side. In the uncoupled representation (|20l) . the operator £^ is diagonal with matrix 
elements given by ^{^ + 1) ^]. The matrix elements of the asymptotic Hamiltonian f^moi 
in basis fl20|) are 

{JMne\{lm,\H^,,\J'M'n'e')\em[) = ( JMHel^^oil J'MW), (22) 

where the expression on the right-hand side is evaluated in Sec. IIA. All that remains to 
complete the definition of the system of coupled equations (12T|) is to evaluate the matrix 
elements of the interaction potential. This is described in the following section. 



dR^ 



+ 2fxE 
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Once the coupled equations are solved, the asymptotic wave function is transformed to 
the field-dressed basis |7)|^m^), which diagonalizes the asymptotic Hamiltonian i^moi- The 
transformation can be written as 



\^)\^m,) = \lm,) CjMnenlJMQe), 



(23) 



JMOe 



where C 



are the components of the eigenvector of Hmoi corresponding to the eigenstate 



7 with energy e^. The matrix of the transformation fl23|) is diagonal in i and m^. The S- 
matrix can be obtained from the transformed wave function using the standard asymptotic 
matching procedure [35|. The cross sections for transitions between the field-dressed states 
of OH can be expressed as 



TT 



a, 



l^r ~ ^2 X^X^ |^77"^«"^m^,m^ ^^imfa' I'm'} 

^ Mtot f.,mi, e',m'f 



(24) 



where the wavevector fc^ = 2fi{E — e^) = 2fj,Eco\\, -Ecou is the collision energy, and the sum- 
mation in Eq. (12^ is performed in a cycle over the total angular momentum projection 

ar constants of OH (in cm'^): Be = 18.55, A = -139.273, 



We used the following molecu 
p = 0.235608, q = -0.03877 |l7L 
(i = 1.68 D was taken from Ref. 



251]. The permanent electric dipole moment of OH(^H) 



17l | . The coupled-channel expansion f[T5]) included all basis 



states with J < y ^ ^ 5, which resulted in a total of 622 coupled channels for Mtot = 2- 
The close-coupled equations (pT!) were solved numerically using the improved log-derivative 



method 



36| on a grid of R from 2 to 65 Qq with a step size of 0.01 Qq. The resulting cross 



sections were converged to within 



70. 



C. Matrix elements of the interaction potential 

The matrix elements of the interaction potential between the states with definite A can 
be expanded in reduced Wigner D-functions 

VAA>{R,e) = J2DoX'-A{o,e,o)Vx,A,-A{R). (25) 

A 

Since in our case |A| = |A'| = 1, only the terms in Eq. ( |25l) with A' — A = 0, ±2 are different 
from zero. They can be obtained by expanding the half-sum and half-difference of the two 
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ground-state potential energy surfaces of A' and A" symmetry 23l. |24| 



-J '^max 

A=0 

-{Va" - Va') = J2 d'o2icos9)Vx2{R) 



(26) 



(27) 



A=2 



where dQ^{cos9) are the reduced Wigner D-functions and Pxij.{cos9) are the associated Leg 
endre polynomials, which are related through 30 1 

"(A-/.)!r/^ 



<(cos^) 



Pam(cos( 



(28) 



.(A + /Z)!. 

The interaction potential (!25|) can be evaluated in the parity-unadapted basis using Eqs. 
(HI), f|T9l) and the generalized spherical harmonics addition theorem [2J, |30| 



c?o/.(cos( 



"^-^(cosl 



2A + 1 



1/2 



A 



«,AO)lAm,(^). (29) 



Combining this expression with Eq. ( l25i) and evaluating the integrals over the products of 
three D-f unctions, we obtain 

( JMfi| (AS I {ime\V{R, 9) \ J'M'Q') | A'S') \i'm'^) = [(2 J + 1) (2 J' + 1) (2£ + 1) {2i' + 1)] 



A,mA 



J A J' 





-m^ mA m'^ 

y<Vx,A'-AiR). (30) 



The expansion coefficients have the property Va,a'-a(-R) = ^a,a-a'(-^) [24]- We note that 
the 3-j symbols in Eq. fl30p vanish unless mx = M' — M = — m'^. Thus, the electro- 
static interaction potential only couples the states with the same total angular momentum 
projection Mtot = M + = M' + m^. A transformation of the interaction potential matrix 
elements ( 1301) to the parity-adapted basis (ITO!) using the symmetry properties of 3-j symbols 
H yields 

{JMne\{^m,\V{R,9)\J'M'n'e')\^'m',) = [(2J + 1)(2J' + 1)(2£ + 1)(2£' + l)]^/2(-)™^+^'-^' 

i \ e 

—rrii mx 






Vx2{R) 



(31) 
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Figure 1: Zeeman energy levels of OH. The initial state for scattering calculations (see text) is 
denoted by the red (light grey) solid line. The values of M for individual magnetic sublevels are 
shown to the right of the curves. 

An analysis of the expression in square brackets shows that the rotational levels in the same 
SO manifold {Vt' = Q) are coupled by the half-sum PES fl26l) . whereas the levels belonging 
to different SO manifolds {Q' = 2 — Q) are coupled by the half-difference PES. Similarly, 
the factor |[1 + ee'{—)^] ensures that the couplings between the states of the same parity 
(e = ±1 ^ e' = ±1) are induced by the anisotropic terms with even A, and those of the 
states of opposite parity (e = ±l^e' = =i=l) are induced by the anisotropic terms with 
odd A. 

III. RESULTS 

The upper panel of Fig. 1 shows the energy levels of OH as functions of the magnetic 
field. At zero field, the absolute ground state of the molecule is a A-doublet with J = 3/2. 
Magnetic fields further split the e and / components of the doublet into four Zeeman levels 
characterized by M = — f? ~|) | (in order of increasing energy). As mentioned above, M 
is rigorously conserved in parallel fields and we will use this label to classify the molecular 
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states. The magnetic levels corresponding to different manifolds cross at i? ~ 0.1 T, where 
the Zeeman shift becomes comparable to the A-doublet splitting. The matrix elements of 
the Zeeman Hamiltonian fllSI) are diagonal in e and independent of its sign. Therefore, the 
e and / Zeeman manifolds are identical and the crossings between the levels of different 
manifolds are not avoided, as illustrated in the upper panel of Fig. 1. Electric fields couple 
the opposite parity states, leading to avoided crossings similar to those encountered in 
molecules 20, 29|. An important difference is that crossings in ^11 molecules occur at 
small magnetic fields S ~ 0.1 T corresponding to typical A-doublet splittings of tenths of 
cm~^. In contrast, the electric-field induced crossings in molecules occur between different 



rotational levels, which requires magnetic fields on the order of several Tesla [2d. l21|. 

The lower panel of Fig. 1 shows the electric field dependence of the e-manifold energy 
levels in the presence of a static magnetic field of 0.01 T. The four e-states shift downwards 
with increasing the field, whereas the corresponding /-states (not shown) shift in the opposite 
direction. It follows from Eq. ( fT6l) that the |M| = | states have larger gf-factors, so their 
energy decreases faster than that of the |M| = | states. As a result, two avoided crossings 
occur at electric fields of about 1 and 1.5 kV/cm shown by the vertical arrows in the lower 
panel of Fig. 1. Note that the location of the crossings depends on magnetic field strength, 
shifting to higher electric fields with increasing magnetic field. In the following, we will 
consider collisions of OH molecules initially in the state \ J = |,M = |,e = —1) denoted by 
the red (light grey) line in Fig. 1. The OH molecules selected in this state gain potential 
energy with increasing magnetic field and can be confined in a permanent magnetic trap. 
The expansion of the initial state in terms of Hund's case (a) basis functions (|23l) is 

|J= |,M = |,e = -1) = 0.985|J= |,n = |,M= |,e = -1) 

+ 0.174|J= |,n = i,M= |,e = -1). (32) 

This equation illustrates that OH is not a pure Hund's ) molecule: different Vt com- 

ponents of the basis ffTOj) are mixed by the cross terms J+S- and J-S+ in Eq. ffTTl) . Because 
the rotational constant of OH is large compared to Zeeman splittings, the mixing coefficients 
in Eq. ( 132|) are independent of M and magnetic field. In the presence of an electric field, the 
initial state (l32l) contains an admixture of basis functions of opposite parity (e = 1), whose 
contribution increases linearly with the field strength. The energy of the state given by Eq. 
( 132|) decreases with increasing electric field as shown in the lower panel of Fig. 1. 
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Collision energy (cm ) 

Figure 2: Cross sections for elastic scattering and spin relaxation in ^He-OH collisions as functions 
of collision energy at different magnetic fields indicated in the graph. The elastic cross section is 
shown for B = 10~^ T. 



Figure 2 shows the cross sections for elastic scattering and inelastic spin relaxation in 
^He-OH collisions as functions of collision energy. The behavior of the cross sections in the 



3l|: the inelastic cross section 



-Ecou limit is dictated by the Wigner threshold laws 

— 1/2 

increases as E^^-^( and the elastic cross section is independent of -EcoU- Magnetic fields 
modify the energy dependence of the cross sections. At small magnetic fields, the cross 
sections continue to decrease with ii^coii down to 10~^ K and start to follow the threshold 
behavior as the energy is further decreased. The turnover point moves to higher energy 
with increasing magnetic field. At very large fields, the spin relaxation cross section always 
increases with decreasing collision energy. This behavior in qualitatively similar to that 
observed for collisions of molecules in S electronic states [11]. The conservation of the total 
angular momentum projection (see Sec. IIB) implies that the spin relaxation transition 



\M 



^) should be accompanied by the transition 



mi + 1 



which leads to a centrifugal barrier in the outgoing collision channel. The centrifugal barrier 
suppresses inelastic processes as long as the ener gy defect between the initial and final 



Zeeman levels does not exceed the barrier height [11|, [3^. At higher magnetic fields (or 
collision energies), the centrifugal barrier is easily surmounted and spin relaxation rates 
increase dramatically, as illustrated in Fig. 2. 
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An interesting feature apparent in Fig. 2 is the rapid increase of spin relaxation as the 
coUision energy is varied through the A-doubhng threshold (0.06 cm~^). This is caused 
by the opening of new relaxation channels in the higher-energy /-manifold (see the upper 
panel of Fig. 1). At collision energies above 0.1 K, both elastic and inelastic cross sections 
display a rich resonance structure. The resonance pattern is rather dense and features both 
shape and Feshbach resonances. This is in contrast with E-state molecules where a single 



20 



29|. The resonances grow in 



(or at most several) shape resonances are typically present 
number with increasing magnetic field. At S = 0.5 T, the spin relaxation cross section shows 
two distinct peaks. We attribute the peaks to shape resonances in the outgoing collision 
channels which are split by the magnetic field. A similar separation of shape resonances has 
been observed for NH(^S^)-He collisions [sll. From Fig. 2, the ratio of the cross sections 
for elastic scattering and spin relaxation varies from 1 to 100 in the temperature interval 
0.01 - IK. Therefore, cryogenic cooling and magnetic trapping of OH using ^He buffer gas 
would be extremely challenging. The rate constant for spin relaxation is 1.2 x 10^^^ cm^/s 
at T = 0.1 K and B = 0.01 T, which corresponds to the OH trapping lifetime of ~0.1 ms 
at the buffer gas density of 10^^ cm~'^. Although spin relaxation is suppressed at collision 
energies below 10 mK and magnetic fields < 0.01 T, this regime is far beyond capability of 
modern cryogenic cooling techniques. 

Equations fl3T|) and fl32l) establish that different Zeeman levels of OH are directly coupled 
by the atom-molecule interaction potential. Therefore, spin relaxation in collisions of ^H 
molecules with ^5*0 atoms is a direct process, in which all Zeeman states get mixed up in a 
collision mediated by electrostatic interactions. The interaction potential for molecules in 
S electronic states is diagonal in spin degrees of freedom, and spin-changing transitions in 

molecules occur via a two-step mechanism through the coupling of the ground and the 
first excited rotational states and the spin-rotation interaction llj. Collision-induced spin 
relaxation in molecules follows a similar mechanism involving the spin-spin interaction. 
In the case of OH, direct couplings of different magnetic sublevels arise due to the anisotropic 
terms (A > 0) in the expansion of the interaction potential over the Wigner D-functions ( l25l) . 

Figure 3 shows the cross sections for spin relaxation as a function of collision energy at 
selected values of the electric field. Electric fields suppress spin relaxation in the ultracold 
regime. The origin of this effect is explained below. At collision energies larger than 0.1 
cm~^, the suppression is much less efficient. Another interesting effect is shifting and splitting 
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Figure 3: Cross sections for elastic scattering and spin relaxation in ^He-OH collisions as functions 
of collision energy calculated for several electric field strengths (in kV/cm) indicated in the graph. 
The magnetic field is fixed at 0.01 T. 

of scattering resonances by electric fields. The suppression of shape resonances is caused 
by the electric field-induced mixing of different partial waves [22]. The number of Feshbach 
resonances grows as the electric field is increased from zero to 10 kV/cm, and the resonances 
shift to lower energies. We attribute this to the electric field-induced couplings between the 
opposite parity states, which are uncoupled in the absence of the field, leading to additional 
avoided crossings and Feshbach resonances. The results shown in Fig. 3 suggest that 
collisional spin relaxation of OH at temperatures below 0.01 K can be suppressed by two 
orders of magnitude with electric fields on the order of 10 kV/cm. The suppression is most 
efficient at low collision energies. The He-OH spin relaxation rate in the absence of an 
electric field is 4.8 x 10~^^ cm^/s at T = 0.01 K and B = 0.01 T. This value decreases 
to 4.4 X 10"^*^ in an electric field of 50 kV/cm. At 0.1 K, the rates are 1.1 x 10"^^ and 
3.4 X 10^^^ cm^/s, respectively. Therefore, cryogenic cooling of magnetically trapped OH 
may be greatly facilitated in the presence of an electric field. A new electromagnetic trap 
for OH molecules 3^ should be particularly suitable for experimental demonstration of the 
electric field-enhanced evaporative cooling. 

In Fig. 4, we show the magnetic field dependence of spin relaxation cross sections at a 
collision energy of 10"^ cm~^ and zero electric field. The cross section shows sharp resonances 
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Figure 4: Magnetic field dependence of spin relaxation cross sections at -EcoU = 10~^ cm~^ in the 
absence of an electric field. 

superimposed on a smoothly varying background. The Feshbach resonances arise due to the 
coupling of the initial channel |J = |,M = |,e = —1) ( l32l) with the closed channels 
I J' = |, M', e = 1) induced by the interaction potential. The closed channels are the bound 
states of the He- ■ ■ OH van der Waals complex that correlate to the upper Zeeman manifold 
in Fig. 1 in the limit R oo. The inset in Fig. 4 demonstrates that resonances may lead 
to a 100-fold enhancement of spin relaxation cross sections at certain magnetic fields. The 
elastic cross section (not shown in Fig. 4) is not affected by Feshbach resonances. 

Figure 5 displays the electric field dependence of spin relaxation cross sections at magnetic 
fields of 10"'^ and 0.01 T and collision energy of 10^^ cm^^. The cross sections increase 
monotonously at low electric fields but exhibit two sudden drops at E > 0.5 kV/cm. The 
location of the dips in Fig. 5 depends on the magnitude of the applied magnetic field. For 
B = 0.01 T, they occur at electric fields of about 1 and 1.5 kV/cm. These values coincide 
with the positions of the level crossings marked by the vertical arrows in Fig. 1. As the 
electric field increases, the energy of the initial state | J = |, M = |, e = —1) becomes lower 
than that of the two \M\ = 1 states. At electric fields larger than ~ 2 kV/cm, both of 
the |M| = i channels become energetically forbidden, and spin relaxation can only occur 
via the |J = |,M = |,e = -1) \ J' = |,M' = -|,e' = -1) transition. An analysis 
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Figure 5: Cross sections for spin relaxation as functions of electric field at i? = 10^'^ T (lower 
curve), B = 0.01 T (upper curve). -EcoU = 1 mK for both curves. The initial state for B = 0.01 T 
is shown by the red (light grey) line in the lower panel of Fig. 1. 

of state-resolved cross sections shows that this transition is the least probable of all spin 
relaxation channels. As a consequence, the total inelastic cross section decreases by four to 
seven orders of magnitude depending on the magnetic field. The data shown in Figs. 2 and 
5 demonstrate that spin relaxation of OH molecules in high electric field-seeking states can 
be completely suppressed by properly chosen combinations of electric and magnetic fields. 
The control is especially robust at low collision energies (on the order of 1 mK) and magnetic 
fields not exceeding 0.01 T. However, Fig. 1 shows that for any given value of magnetic field, 
it should be possible to chose an electric field at which the |M| = | channels are closed. The 
necessary electric field can be determined from the positions of the crossings in the lower 
panel of Fig. 1. 

IV. SUMMARY 

We have developed a rigorous quantum theory for collisions of molecules in ^H electronic 
states with structureless atoms in the presence of superimposed electric and magnetic fields. 
The matrix elements of the electrostatic potential and molecule-field interactions have been 
derived in the fully uncoupled representation of Hund's ) basis functions. The theory 
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has been applied to elucidate the mechanisms of inelastic transitions in low-energy collisions 
of OH molecules with He atoms. Our results suggest that spin relaxation in collisions 
of OH molecules proceeds via direct coupling of different Zeeman states induced by the 
anisotropy of the atom-molecule interaction potential. The rate constants for spin relaxation 
at temperatures above 0.1 K are on the order of 10^^^ cm^/s, leading to very short trapping 
lifetimes. We conclude that sympathetic cooling of OH molecules using cryogenic '^He gas 
at densities > 10^^ cm~^ and temperatures 0.1 -IK appears unfeasible. 

We have demonstrated that spin relaxation of OH molecules at temperatures below 0.01 K 
can be efficiently manipulated by electric fields of moderate strength (~10 kV/cm) available 
in the laboratory. The mechanism of control is based on suppressing certain relaxation 
pathways using superimposed electric and magnetic fields. This technique can be used to 
facilitate evaporative cooling of ^H molecules in electromagnetic traps (34|. We have found 
that electric fields modify the collision energy dependence of inelastic cross sections and 
may lead to the formation and splitting of Feshbach resonances. The magnetically tunable 
Feshbach resonances may be used to create weakly bound He- ■ ■ OH complexes. It would be 
interesting to explore the effects of magnetic fields on spin-orbit, vibrational, and rotational 
predissociation of these complexes. The methods to control inelastic collisionspresented 
in this work may be realized experimentally using Stark decelerated beams |5|, [l2| and 
electromagnetic traps 0, Q. 
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